Pyrochlore Photons: The U(l) Spin Liquid in a S = 1/2 Three-Dimensional Frustrated 
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We study the S = 1/2 Heisenberg antiferromagnet on the pyrochlore lattice in the limit of strong 
easy-axis exchange anisotropy. We find, using only standard techniques of degenerate perturbation 
theory, that the model has a (7(1) gauge symmetry generated by certain local rotations about the z- 
axis in spin space. Upon addition of an extra local interaction in this and a related model with spins 
on a three-dimensional network of corner-sharing octahedra, we can write down the exact ground 
state wavefunction with no further approximations. Using the properties of the soluble point we 
show that these models enter the U(l) spin liquid phase, a novel fractionalized spin liquid with an 
emergent U(l) gauge structure. This phase supports gapped S z — 1/2 spinons carrying the U(l) 
"electric" gauge charge, a gapped topological point defect or "magnetic" monopole, and a gapless 
"photon," which in spin language is a gapless, linearly dispersing S z = collective mode. There are 
power-law spin correlations with a nontrivial angular dependence, as well as novel U(l) topological 
order. This state is stable to all zero-temperature perturbations and exists over a finite extent of 
the phase diagram. Using a convenient lattice version of electric-magnetic duality, we develop the 
effective description of the U(l) spin liquid and the adjacent soluble point in terms of Gaussian 
quantum electrodynamics and calculate a few of the universal properties. The resulting picture is 
confirmed by our numerical analysis of the soluble point wavefunction. Finally, we briefly discuss 
the prospects for understanding this physics in a wider range of models and for making contact with 
experiments. 



I. INTRODUCTION 

The search for quantum spin liquid states in frus- 
trated magnets can be traced back at least as far as 
the early suggestion of a resonating valence bond state 
in the triangular lattice Heisenberg model 1 . Almost 15 
years later, the suggestion of such a state in the un- 
doped high-T c cuprates^ set off an explosion of inter- 
est in two-dimensional spin liquids (i.e. Mott insulators 
at half- filling with no broken symmetries). Frustrated 
Heisenberg models on the square, triangular and kagome 
lattices have all received significant attention as can- 
didate systems for quantum disordered ground states. 
While there has been comparatively little theoretical 
work on quantum spin liquids in three-dimensional frus- 
trated magnets, materials with magnetic ions on the py- 
rochlore lattice (Fig. [2| may be good candidates for spin 
liquids and other exotic states. To give one example, 
recent neutron scattering experiments on ZnCr204, a 
S = 3/2 pyrochlore Heisenberg antiferromagnet, suggest 
a nontrivial disordered stated above a transition to Neel 
order accompanied by a lattice distortion! at 12.5 K. 

Meanwhile, much work has been devoted to under- 
standing the properties of possible spin liquid states, 
independent of their existence in particular microscopic 
models. Most of the proposed spin liquid states support 
deconfincd S = 1/2 spinons; such states are fraction- 
alized in that some of the elementary excitations carry 
quantum numbers that are fractions of those allowed in a 
finite-size system. Fractionalized states can be precisely 
characterized by their topological orders, which in the 
simplest scenario is associated with the topological sec- 



tors of an emergent, deconfining Z2 gauge field 6 . The 
"vison," a gapped vortex-like excitation that carries the 
Z2 flu3#i, must also be present. While topological order 
does not require a liquid ground state and can coexist 
with conventional long-range order, we believe it is prob- 
ably common in spin liquids and hence in some of the 
nearby ordered states. Very recently, many of these ideas 
have been put on firmer ground by the emergence of sev- 
eral microscopic models supporting stable fractionalized 
phases in two and three dimensions&2*i2iii*i2iiiLii. 

Despite these recent theoretical successes, an unam- 
biguous experimental realization of these ideas is still 
lacking. Indeed, spin liquid states seem rather rare; 
is topological order rare as well? For ^-fractionalized 
states this question is difficult to answer, because the 
gapped visons have no effect on easily measurable low- 
energy properties. Clever proposals have been madei^ 
and carried ou t 16 i 17 to directly detect topological order 
in the cuprates (with negative results thus far), but these 
experiments are difficult and rely on properties of the 
phases proximate to a topologically ordered state. Z2 
topological order is difficult enough to observe that it is 
impossible to say at present how rare or common it is. 

Fortunately it may be possible to shed new light on 
the experimental situation. In this paper, we present two 
models of three-dimensional S — 1/2 frustrated magnets, 
one on the pyrochlore lattice, the other on a related net- 
work of corner-sharing octahedra (the links of the cubic 
lattice, as shown in Fig. EJ- Both of these models exhibit 
a novel fractionalized phase, the U(l) spin liquid. This 
state has an emergent U(l) gauge structure that gives 
rise to several remarkable properties: there is a gapless 
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"artificial photon" excitation, a gapped spinon carrying 
"electric" gauge charge, a gapped "magnetic" monopole, 
an emergent 1/r "Coulomb" potential between pairs of 
spinons and monopoles, and novel £7(1) topological order. 
If this phase exists in a real material, the gapless photon 
should have important implications for low-energy ther- 
modynamics, transport and spectroscopy; therefore £7(1) 
fractionalization may be easier to find in experiments. 
Such states, thus far realized in large- N spin models^ 
and bosonic Hubbard-type modelsiSiiiiiS, arise as the de- 
confined or Coulomb phase of compact £7(1) lattice gauge 
theory. While most work on spin liquids has focused on 
d < 2, motivated by the cuprates and the conventional 
wisdom that quantum fluctuations are more effective at 
destroying long range order in low dimensions, the £7(1) 
spin liquid only occurs in d > 3; for d < 2 the Coulomb 
phase of compact £7(1) gauge theory (with gapped mat- 
ter) is always unstable due to instanton effects^. 

Both models are of intrinsic interest as examples of 
tractable but nontrivial frustrated magnets. The py- 
rochlore model is particularly appealing due to its sim- 
plicity: its derivation begins with the nearest-neighbor 
5 = 1/2 Heisenberg antiferromagnet. Taking the limit of 
large easy-axis exchange anisotropy J z Jj_ simplifies 
the problem by breaking the spectrum into extensively 
degenerate manifolds with large separations of 0{J Z ). It 
is then possible to write an effective Hamiltonian describ- 
ing the splitting of the low-energy manifold, using stan- 
dard techniques of degenerate perturbation theory in J± . 
This effective Hamiltonian has a £7(1) gauge structure, 
which forms the foundation for our subsequent analysis 20 . 
Another point of view, equivalent at the level of pertur- 
bation theory but perhaps with broader implications in 
more general scenarios, is that the low-energy sector of 
the model is unitarily equivalent to a £7(1) gauge theory. 
It is not obvious how to treat the resulting model analyt- 
ically, but upon addition of an extra six-site interaction 
term it can be tuned to a soluble point where it is pos- 
sible to write an exact ground state wavefunction with 
no further approximations. The models can be reinter- 
preted as quantum dimer models (QDMs), and the extra 
term as the analog of the Rokhsar-Kivelson (RK) poten- 
tial in the square lattice QDM 21 . As will be explained in 
detail below, the properties of the soluble point allow us 
to locate the U(l) spin liquid adjacent to it. Since this 
state is stable to all zero-temperature perturbations, it 
persists over a finite extent of the phase diagram (Fig-QJ. 
Furthermore, stability to large but finite J z implies that 
the £7(1) gauge structure persists in the absence of mi- 
croscopic local symmetries and is truly emergent. On 
the purely theoretical side, we believe these models give 
the first examples of £7(1) gauge theories that have a de- 
confining phase even in the limit of infinitely strong bare 
coupling. The first such Zi gauge theory was discovered 
only recently by Moessner and Sondhi&. 

The effective theory of the £7(1) spin liquid and the 
soluble RK point is simply Gaussian quantum electrody- 
namics (QED). At the RK point, which is itself a spe- 
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FIG. 1: Phase diagram for both models. The parameter 
V/Jring is the relative strength of the Rokhsar-Kivelson po- 
tential and the XY ring exchange that obtains in the easy-axis 
limit of the Heisenberg model. The soluble point is located 
at V/Jrin S = 1, which is a special deconfined point of the ad- 
jacent 17(1) spin liquid. Just to the right of the soluble point 
the models go into an Ising ordered state. Sufficiently far to 
the left we expect Ising order, while at intermediate values 
of V/Jring states with broken translation symmetry but no 
magnetic order are also possible. Immediately to the left of 
the soluble point, the £7(1) spin liquid exists over a finite (but 
unknown) extent of the phase diagram. 



cial deconfined limit of the generic phase, the "electric 
stiffness," or coefficient of E 2 in the Hamiltonian, van- 
ishes. This is a higher-dimensional generalization of the 
effective picture of the square lattice QDM in terms of a 
coarse-grained height field^. 

The 17(1) spin liquid has power-law correlations with 
nontrivial angular dependence, novel U (1) topological or- 
der, and supports gapped S z = 1/2 spinons, a gapped 
topological point defect (the "magnetic" monopole) and a 
gapless S z = collective mode corresponding to the pho- 
ton of the gauge theory. The latter excitation makes an 
additive T 3 contribution to the low-temperature specific 
heat, and should affect various other low-energy proper- 
ties of £/(l)-fractionalized phases (either the £7(1) spin 
liquid, or phases with coexisting conventional and topo- 
logical order). If such a phase exists in a real material, 
we speculate that it may be possible to probe "photons" 
with photons via Raman scattering. 



A. Outline 

We begin Sec. [H] with a derivation of the pyrochlore 
model starting from the Heisenberg antiferromagnet. In 
Sec. Ill Bl the cubic (or corner-sharing octahedra) model 
is discussed. The remainder of Sec. [H] is concerned with 
demonstrating the equivalence of the spin models to frus- 
trated compact 17(1) gauge theories, and developing a 
useful lattice version of electric-magnetic duality. 

Beginning from the dual description, Sec. IIHI develops 
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FIG. 2: The pyrochlore lattice (left), and one up-pointing 
tetrahedron (right). One sublattice of tetrahedra is shaded, 
and the other transparent. The thickened bonds show the 
location of a pyrochlore hexagon. Each such hexagon is a 
member of one of four orientations of kagome lattice planes. 
The numbering of sites in the up-pointing tetrahedron on the 
right is the convention used in the text. For i = 0, 1,2, the 
fee Bravais lattice vector &i points in the direction given by 
looking from site 3 to site i. 
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FIG. 3: Depiction of the processes contributing to the third- 
order degenerate perturbation theory for the easy-axis py- 
rochlore Heisenberg antiferromagnet. Processes (A) and (B) 
give only trivial constant shifts of the energy. Process (C) 
leads to an XY ring exchange term acting on hexagonal pla- 
quettes. 



the effective description of the U(l) spin liquid and the 
soluble point in terms of Gaussian quantum electrody- 
namics. Corrections to effective action and to the scaling 
equalities between microscopic and effective degrees of 
freedom are discussed in Sec. IIII CI Sec. IIVI contains a 
discussion of the universal properties of the U(l) spin liq- 
uid, including its novel U(l) topological order. In Sec.lVI 
we present our analysis of the soluble point ground state 
wavefunction, which gives strong support for the valid- 
ity of our effective picture. We conclude in Sec. IVII with 
a discussion of open issues, focusing on the challenging 
problems of understanding this physics in a broader range 
of models and looking for £/(l)-fractionalizcd phases in 
real materials. 



II. MODELS AND MAPPINGS 



A. Pyrochlore Model 

We begin with the nearest-neighbor S — 1/2 Heisen- 
berg antiferromagnet on the pyrochlore lattice. This 
structure is a three-dimensional network of corner- 
sharing tetrahedra (Fig. |2J). It can be obtained by trans- 
lating one "up-pointing" tetrahedron (shown on the right 
of Fig. |5J) through the fee Bravais lattice vectors R = 
noSLo+niai+n2a 2 . We choose a = x, ai = x/2+V3y/2, 
and a 2 = x/2 + y/2^/3 + ^/2/3z. Basis vectors for the 
reciprocal lattice are defined by = y/2~Treijka.j x a^, 
so that &i ■ hj = 2TrSij. The four sites in each unit cell 
are distinguished by an index i = 0, . . . , 3, as indicated 
in Fig. [3 Lattice sites are denoted either by single italic 



letters like i, or by pairs (BL, i) when we wish to specify 
the position of a site within the unit cell. 

Up to a constant the Hamiltonian can be written as a 
sum over tetrahedra: 



(i) 



where St = X)iet * s ^ ne total spin on the tetrahedron 
t. Following the analysis of a generalized kagome Heisen- 
berg antiferromagnet in Ref. [9j, we introduce easy-axis 
exchange anisotropy: 



H = 
Hi = 

w = 



Ti.i + TL' , 



(2) 
(3) 

(4) 



(ij) 



where J z 3> J±. This reduces the global SU(2) invari- 
ance to U(l) x We first consider the point J± = 0, 
where H reduces to a classical Ising model, with ground 
states specified by S* — on all tetrahedra. It was ar- 
gued by Anderson^ that, almost identically to Pauling's 
model for water ice24, this Ising model has an extensive 
ground state degeneracy (i.e. finite T = entropy per 
site). 

A small J± > introduces quantum fluctuations and 
lifts the extensive degeneracy; this splitting is encapsu- 
lated in an effective Hamiltonian using standard tech- 
niques of perturbation theory. The first-order contribu- 
tion is easily seen to vanish. We will need to go to third 
order, where we have the general expression: 
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V V V 

H'—H'+H'—H'—H' 

rii rti rii 



(l-V). (5) 
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Here V projects onto the orthogonal complement of the 
ground state manifold. To describe the processes con- 
tributing in Eq. 0, it is useful to work in the standard 
hardcore boson language for the spins, where S z — ±1/2 
corresponds to the presence/absence of a boson. Each 
term in TL' hops bosons along nearest-neighbor bonds; 
acting on a state in the low-energy manifold, each hop 
creates two tetrahedra with S% ^ 0. At second order 
in TL' , bosons can hop and then return along the same 
bond (Fig. [3JY) . This can always occur on 4 bonds in ev- 
ery tetrahedron, thus giving only a constant contribution 
to the energy. At third order another constant contribu- 
tion arises from single bosons (or holes) hopping around 
triangular faces (Fig. 03) . There is also a nontrivial ring 
exchange process acting on the hexagonal plaquettes (see 
Fig. I2J , where hexagons containing three evenly spaced 
bosons can be rotated as shown in Figure|3p. The result- 
ing effective Hamiltonian is 

H efJ = (Jl/J z )(J ± /J z -l)N t (6) 

+ Jring'^jiS^ S 2 S3 S 4 S 6 +h.C.), 

O 

where N t is the total number of tetrahedra, J T ing = 
3J^/2J^, and the sum is over hexagonal plaquettes. The 
labelling of the spin operators inside the sum is given by 
moving around each hexagon in an arbitrary direction. 
Note that \H e ff,Sf] — 0, as must be true for any ef- 
fective Hamiltonian acting in the low-energy manifold, 
whatever the form of TL' . 

We focus on the extreme easy axis limit described by 
TL e f f , but note in passing that a finite but large J z would 
introduce small fluctuations out of the ground state man- 
ifold. While these will not affect universal properties, 
they can matter for the short-distance correlation func- 
tions of some microscopic operators. This can be under- 
stood formally by a more sophisticated execution of the 
perturbation theory in J± that accounts for splitting of 
the low-energy manifold and mixing of higher states on 
an equal footing 2 ^. The main result is that the problem 
at finite J z can be mapped, by a unitary transforma- 
tion, order-by-order in J± onto a transformed Hamilto- 
nian acting only within the low-energy manifold where 
— 0. This mapping accounts for finite J z by gener- 
ating nontrivial perturbative relations between physical 
and transformed operators; however, we ignore these cor- 
rections for simplicity and use only the results of the stan- 
dard degenerate perturbation theory described above. 

It is possible, and will be convenient, to change the sign 
of the ring term by a similarity transformation. On any 
given site we can make the transformation S z — ► S z and 
S — > — by making a 7r-rotation about the z-axis in 
spin space. One transformation with the desired effect, 
consisting of 7r-rotations on a pattern of sites, is: 

Si - S? (7) 
S± - exp^Q, • R)S± , (8) 

where Qo = Qi = (bi + b 2 )/2 and Q 2 = Q3 = 0. 




FIG. 4: A small piece of the diamond lattice. The links form 
hexagonal loops corresponding to the pyrochlore hexagons. 
These are the shortest possible closed paths on the diamond 
lattice. The hexagon with three thickened bonds depicts the 
dimer positions on a flippable hexagon. The alternating full 
and empty bonds correspond to alternating up and down 
spins. 



After this transformation the Hamiltonian takes the 
form 

TL p = — J ring {S^ S 2 S3 S4 S^ Sg + h.C.), (9) 

o 

where the constant terms have been dropped. Mod- 
els similar to this one on the kagome^, squar o 2 9i 2 7 , 
triangular— other lattices^, where XY ring exchange of 
spins or bosons is a dominant term, have recently been 
shown to exhibit a variety of unusual phases and crit- 
ical behavior. The physics of the pyrochlore ring ex- 
change model should be accessible to quantum Monte 
Carlo studies; while the original Hamiltonian in Eq. @ 
has a sign problem, Ti v does not. 

Tip can be reinterpreted as a quantum dimer model 
on the diamond lattice (Fig. Q}, with two dimers touch- 
ing every site. To see this, observe that the centers of 
the pyrochlore tetrahedra form a diamond lattice. Each 
nearest-neighbor diamond link passes through exactly 
one pyrochlore site, so we can reinterpret the pyrochlore 
spins as diamond link variables. The smallest closed 
loops in this lattice contain six links and correspond to 
the pyrochlore hexagons. We say a dimer is present on a 
given bond if S? = 1/2, or absent if Sf = -1/2. = 
becomes the constraint that every diamond site touches 
two dimers, and the ring exchange move is the most lo- 
cal dynamics preserving this constraint. Each term in TLp 
acts on a "flippable" hexagon, one containing alternating 
full and empty bonds as in Fig. by rotating the dimers 
around it. Non-flippable hexagons are annihilated. 

As first realized by Rokhsar and Kivelson, dimer mod- 
els generically have a point in their parameter space 
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where it is possible to write down the exact ground state 
wavefunctior. 21 . To reach this point in our model, we 
add the term Ti.v — VNf, where Nf is the number of 
flippable hexagons. The RK point obtains for T-Lrk = 
Hp + JringNf (i.e. V = Jring), and the ground state is 
an equal-weight superposition of all possible dimer cov- 
erings of the lattice that satisfy the constraint of two 
dimers touching every site. In the spin language, this 
wavefunction can be written as the projection of a trans- 
verse ferromagnet: 

\1>RK) = {l-V)l[\Sf = 1/2), (10) 

i 

where, as in Eq. JSJ, (1 — V) projects onto the — 
manifold. For completeness, we also express Nf in terms 
of spin operators: Nf = J2o P fiw(°)' where Pfu p (o) 
gives unity acting on a flippable hexagon and zero other- 
wise. One has: 

CT =±i jeo 

We will be interested in the properties of the generalized 
ring model H. p + TLv in the vicinity of the soluble point. 

B. Cubic Model 

Largely to simplify the geometry of the presentation, 
we introduce an alternate model that we find has many of 
the same properties as its pyrochlore analog. The model 
is the QDM on the cubic lattice, with three dimers touch- 
ing every site. We consider only the most local dynam- 
ics, which rotates the configuration on square plaquettes 
with two dimers on opposite sides, and the corresponding 
Rokhsar-Kivelson potential that counts flippable squares. 
Reversing the mapping above, we can also think of this 
as a spin model with S = 1/2 on the links of a cubic lat- 
tice, or, equivalently, on the sites of a lattice of corner- 
sharing octahedra with their centers at the cubic sites 
(Fig. [5J . The octahedra play the role of the pyrochlore 
tetrahedra, with the total spin on each S* ct — 0. We 
denote cubic sites by boldface letters like r and identify 
the links by specifying either pairs of adjacent sites, or 
one site and the direction of the link. For example, the 
link connecting a site r with its nearest neighbor in the 
positive x-direction (r + x) is denoted by (r, r + x) or 
(r,x). Using this notation, we express the dimer kinetic 
term as a 4-site XY ring exchange for the spins: 

n c = -j rmg J2(sts^s+s^ + h. c .) (12) 

□ 

r 

where the numbering in the first line runs around the 
perimeter of each square plaquette, and in the second 
line only one orientation of plaquette is shown explicitly. 




FIG. 5: Illustration that the links of the cubic lattice are 
equivalent to the sites of a lattice of corner-sharing octahedra. 



As before, we will be interested in the vicinity of the 
soluble point of H c + Hv 

It is also interesting to note that, as for H. p , the pure 
ring exchange model 7i c can be derived as the easy-axis 
limit of a Heisenberg antiferromagnet. In this case one 
begins with spins on the network of corner-sharing octa- 
hedra with nearest-neighbor exchange, and an additional 
exchange of the same sign and magnitude between spins 
at opposite points of the same octahedron. The Hamil- 
tonian can be written as a sum over octahedra and the 
analysis proceeds exactly as before. 



C. Frustrated U(l) Gauge Theory 

Both the cubic and pyrochlore models have an ex- 
act U(l) gauge invariance, as is generally the case in 
dimer models. In the pyrochlore model, this arises be- 
cause of the local integer conserved quantity The 
local symmetry is generated by rotations about the z- 
axis in spin space on all the sites in a given tetrahedron: 
G t (a) = exp(iaSt). These generators commute with one 
another and the Hamiltonian. It is important to empha- 
size that this gauge symmetry is generated by physical 
transformations and has nothing to do with any redun- 
dancy in our description. These statements also hold for 
the cubic model, with octahedra substituted for tetrahe- 
dra. For now, we shall focus on the cubic model for ease 
of presentation. All of the statements in this and the 
following section can be generalized simply to the py- 
rochlore case; this is summarized in Section 111 El below. 

The gauge structure suggests we may gain insight by 
thinking about the model as a U(l) lattice gauge theory, 
an approach that proved helpful in understanding the 
square lattice QDM 2 ^. In fact, it is possible to formally 
rewrite TL C as a pure gauge theory on the cubic lattice. 
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To see this, we find it convenient to soften the hardcore 
constraint on the bosons by explicitly introducing an on- 
site repulsion term in the Hamiltonian. We also go to 
quantum rotor variables rt rr < G Z and rr ' G [— 7r, 7t), 
with commutation relations [4> r i,n r ij] — iSijS rr i. (Here 
and elsewhere we commit a standard abuse of notation 
and formally denote a constraint on the eigenvalues of 
an operator as a constraint on the operator itself.) Using 
the relations S* ~ n — 1/2 and = exp(±i0), we have 

«c = ^5> rr ,-l/2) 2 (13) 

<rr'> 

- JT^COs(0i - 2 + (f> 3 - 4 ), 
□ 

where the numbering inside the cosine proceeds around 
the perimeter of the given square. This is a faithful 
representation of the spin model Eq. (|12|) in the limit 
U/K — > oo, which just imposes the hardcore constraint 
n rr > = 0, 1. 

We now define an orientation on the cubic links, which 
we take to point out of the A sublattice and into the B 
sublattice. This allows us to define oriented link variables 
by: 

e rr > = ±(n rr , - 1/2) (14) 

= ±0rr'- (15) 

Here we take the plus/minus sign when r lies in the A/B 
sublattice. Since e rr ' = — e r ' r , these variables can be 
thought of as components of vector fields taken along 
the links of the lattice. Putting these definitions into the 
Hamiltonian, we have 

(rr'> □ rr'GD 

The sum inside the cosine is taken in an oriented fashion 
around the links of the given square plaquette, and is thus 
a discrete line integral of a rr / along a smallest possible 
closed path. This is a lattice version of the curl, so we 
define: 

(curl a) □ = ^a rr /. (17) 

rr'eD 

Eq. (|16J) is invariant under gauge transformations writ- 
ten in the usual form a rr > —>■ a rr > +Xr' — Xr- Therefore a rr ' 
behaves like a vector potential, and since [a rT >,e Tr /] = i, 
e rr > plays the role of an electric field. Since the vector 
potential is a 27r-periodic phase this is evidently a com- 
pact gauge theory. The "electric" charge is given by the 
lattice divergence of e rr > : 

(dbre) r =J2 e rr> = ±S* oct , (18) 

r'< — r 

where the sum is over the sites adjacent to r. In the 
ground state there is no gauge charge since S% ct = 0, and 
single gauge charges have a large gap of order J z . 



Suppose r is a site in the A sublattice. Acting with S+ x 
creates two octahedra with S% ct — 1, at r and r+x. In the 
gauge theory language, there is a positive electric charge 
at r and a negative one at r + x, where the sign of the 
charge comes from the orientation convention Eq. (T3}. 
We can now act with the infinite string operator: 

oo 

Ostring = \\ ^[r+(ra+2)x+ny],3/^[r+(ra+l)x+ny],a; - (^) 
n=0 

This hops the gauge charge originally at r + x off to infin- 
ity, leaving an isolated S z = 1 octahedron at r. The two 
S* = 1 octahedra together carry S^ otal = 1 , so the single 
remaining octahedron is evidently a S z — 1/2 spinonM^. 
The spinons are single electric gauge charges and can 
propagate freely in a deconfining phase of the gauge the- 
ory; such a state is therefore fractionalized. We note 
that, because spinons cannot hop from one sublattice to 
the other, there are in fact two spinon flavors. 

Because the theory is compact, magnetic charge is 
also allowed. Define dual lattice sites by serif charac- 
ters r = r + (x + y + z)/2; these are located at the centers 
of the cubic "boxes" of the direct lattice. The links of 
the dual lattice are naturally associated with the square 
plaquettes of the direct lattice. Using this correspon- 
dence we define a magnetic field on the dual links by 
7r6 rr / = (curla)p. The sense of the lattice curl is taken 
counterclockwise looking from the point r' to r. The lat- 
tice divergence (div6) r gives the magnetic charge inside 
the box at r. Naively this divergence vanishes since each 
term a rr ' occurs twice, with opposite signs, but this is 
not the case because 6 is a periodic variable invariant un- 
der 6-^6 + 2. In fact we have (div6) r = 2n r for integer 
n r ; the magnetic charge is automatically quantized. It is 
convenient to take 6 rr < G [—1, 1), so that n r = 0, ±1,±2 
on each box. It is also possible to have n r = — 3 in our 
convention, but this measure zero point in the configu- 
ration space of magnetic fields should be ignored. While 
electric charges are locally conserved in the low-energy 
manifold, there is no such conservation law for magnetic 
charge and we expect n r ^ even in the ground state. 
While this means the ground state always has some local 
fluctuations of magnetic charge, it does not necessarily 
contain monopoles, which are smoothly- varying defect 
configurations unaffected by a small amount of coarse- 
graining. 

The model Eq. I|16|) looks identical to the standard 
Hamiltonian formulation of compact U(l) lattice gauge 
theory, but there is one difference of critical importance 
in the limit of interest U/K — > oo. Here, the electric field 
takes on half-integer rather than integer values. In the 
case of integer electric fields, the vacuum in the large-[7 
limit is trivial: e — everywhere with small fluctuations. 
In our model this limit enforces the nontrivial constraint 
e = ±1/2. This is another expression of the inherent 
frustration, so we refer to the model as a "frustrated 
gauge theory." This situation is essentially the same 
as that arising in the large-iV limit of bipartite SU(N) 
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antiferromagnets2i*2£, and in the gauge theory descrip- 
tion of the square lattice QDM 29 . There, one considers 
an integer-valued electric field in a background of static 
charge; of course, our model can also be viewed this way 
by making a shift of the electric field. 

To gain insight into the possible phases of TL C , we 
briefly review the properties of the standard unfrustrated 
gauge theory, which has the same Hamiltonian as our 
model but an integer electric fieldS. There are two 
phases, separated by a transition at (U/K) c » 1. For 
U/K > (U/K) c , the "strong coupling" side of the tran- 
sition, the model enters a confining phase smoothly con- 
nected to the trivial U = oo vacuum e rr ' = 0. In this 
phase all excitations are gapped, and static sources of 
the electric field are confined by a linear potential. For 
U/K < (U/K) c one enters the deconfining Coulomb 
phase, so named because static gauge charges interact 
via a 1/r Coulomb potential and are thus free to prop- 
agate. At low energies, the effective description of the 
Coulomb phase is simply Gaussian QED, so there is a 
gapless, linearly dispersing photon with two transverse 
polarizations. 

In the Coulomb phase, magnetic monopoles are gapped 
and interact via a 1/r magnetic Coulomb potential. The 
monopoles incorporate the compact nature of the mag- 
netic field, which is not important at low energies in the 
deconfined phase. In the confined phase, however, the 
magnetic field fluctuates wildly, its periodicity is impor- 
tant, and the monopoles have proliferated and condensed. 
This distinguishes the two phases: in the Coulomb phase 
the monopole propagator decays exponentially, while in 
the confined phase it goes to a constant. This distinction 
is more robust than the Wilson loop, which fails to dif- 
ferentiate between the phases in the presence of matter 
fields. 

Returning to the frustrated gauge theory, we can infer 
that it also has a Coulomb phase at small U/K. This 
should be so because in the deconfined phase the discrete 
nature of the electric field is unimportant, so the half- 
integer nature of e rr ' will not play a role. We will be 
interested in whether the Coulomb phase survives in the 
opposite limit of strong coupling, perhaps stabilized by 
additional terms in the Hamiltonian. 



D. Electric-Magnetic Duality 

In recent work on other models with XY ring exchange, 
it has been useful to make a novel "plaquette duality" 
transformation 26 28 . In the case of the frustrated gauge 
theory, plaquette duality is in fact identical to the more 
familiar electric-magnetic duality for U(l) gauge theories. 
We shall make use of the dual cubic lattice defined in the 
preceeding section. We define oriented variables on the 
dual links a„i € 7rZ and b rr > £ [—1,1), which we take to 
be canonically conjugate: [6, a] — i on the same dual link, 
zero otherwise. As discussed above, b rr ' is the magnetic 



field, defined by: 

irb rr > — (curla)o, (20) 

where the sense of the curl is again taken counterclock- 
wise looking from r' to r. The conjugate variable will 
play the role of a "dual vector potential," and is related 
to the electric field by: 

7r(e rr / -e° r ,) = (curl <*)□., (21) 

where the curl is taken around the dual square plaquette 
encircling the direct link joining r to r', and again the 
orientation is given by the right hand rule. Here e® r , = 
±1/2 is a static, divergenceless background field; this is 
necessary for consistency with e rr / € Z+l/2. With these 
definitions, the dual commutation relations are consistent 
with the original ones. 

In dual variables the Hamiltonian is 

Wed = ^X (( curla b* +ne° r ,) 2 - if ^cos(vr& rr /), 

<rr') 

where the first sum is over dual plaquettes. The con- 
straint of zero electric charge is now automatically sat- 
isfied, since dive = div(e° + curla/7r) = 0. How- 
ever, magnetic charge can take on continuous values 
in the dual variables, and one should impose the con- 
straint div b £ 2Z. This constraint commutes with TL c d 
since it is invariant under the dual gauge transformations 
a„> — > a„> + A r < — A r , where A r € 7rZ. Note that because 
magnetic charge is not locally conserved, we are not al- 
lowed to demand div b = 0. 

It is useful, and enlightening, to write down the Eu- 
clidean action obtained by a Trotter expansion in eigen- 
states of the dual vector potential. As usual, one begins 
with the partition function 

Z cd = Tr(exp(-/?H C(i )p). (23) 

Here V — Y[ r V r is a projection operator imposing the 
quantization of magnetic charge: 

V r = <K(div&) r -2a rr /7r) ( 24 ) 

= 2 XI exp (* a rr(div6) r ). 

Breaking the exponential in Eq. lll'.'il) into N T = f3/e 
timeslices, and inserting V once in every timeslice, one 
has: 

Zed = Tr(exp(-ett c d)P) (25) 
= e MScd), (26) 

where the index fj, = t,x, y, z. The imaginary time com- 
ponent of a ^ comes from the Poisson-resummed form of 
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the projector Eq. J33J. The spatial components enter 
as the eigenvalues of the states |{a r i(r)}) used to form 
the resolution of the identity at each timeslice. Follow- 
ing very similar manipulations to those in Appendix A 
of Ref . I|26|) , one obtains the dual action: 

S cd = ^ln(J^) J2 (Arttrr' + arr - «rv) 
T,(rr') 

+ Y~2 E (W lc 0n* +™° rr ,y, (27) 

n r,P* 

where A T / = /(r + e) - /(r). 

The action Eq. I|27(l is essentially a higher-dimensional 
generalization of the height model partition function ar- 
rived at by similar manipulations in the context of the 
square lattice QDM2&. Significantly, it differs in having 
a local rather than a global invariance, under spacetime- 
dependent dual gauge transformations of the form a M — > 
ct M + A^A. In fact, S c d has the same structure as non- 
compact lattice QED, except for the discrete nature of 
the fields. This encodes the important physics of the 
magnetic monopoles. In Sec. IIIII below, it will be use- 
ful to imagine softening the constraint of discreteness on 
the fields to interpolate between the dual partition func- 
tion and an effective description of the Coulomb phase. 
With the soft constraint (implemented by the "correc- 
tions" in Eq. (|55Jl ). the theory is identical to a more 
familiar dual representation of U(l) gauge theory con- 
sisting of a non-compact gauge field minimally coupled 
to scalar monopoles. 



then takes the form of the diamond lattice frustrated 
gauge theory: 

U P = | E e rr< - K E C0S ( E ■ ( 29 ) 

<rr') O rr'eO 

It is evident that the lattice curl now naturally lives on 
the hexagons of the diamond lattice. 

Again the electric charge has a simple interpretation 
in the spin language: 

(dive) r = ±Sf. (30) 

Tetrahedra with S* = ±1 are now the S z = 1/2 spinons 
carrying unit gauge charge. Single spinons can be created 
by a string operator similar to Eq. (|19fl . For small U/K 
the model should again enter a deconfining phase where 
the spinons are free to propagate. 

We define a dual lattice of plaquette variables by 
putting a site at the center of every pyrochlore hexagon. 
This is also a pyrochlore lattice, and it will be useful to 
think of its sites as the links of a dual diamond lattice 
with sites labelled by serif characters r. Each hexagon of 
the dual lattice encircles a link of the direct lattice, and 
vice versa. As before, magnetic charge lives on the dual 
lattice sites. The dual variables, again with the commu- 
tator [b, a] = i on the same link, are defined on the dual 
links by: 

7r(e rr ' — e® r ,) — (curla)o* (31) 
irb n ' — (curla)o, (32) 



E. Pyrochlore Gauge Theory and Duality 

We now return to the pyrochlore ring model TL p . In 
this case the diamond lattice with sites at the centers of 
the tetrahedra (discussed in Sec. Ill A|) plays the role the 
cubic lattice did for the cubic model. Denoting diamond 
sites by boldface characters, we soften the hardcore con- 
straint on the bosons and go to quantum rotor variables 
living on the diamond links: 

K P = |]>> rr ,-l/2) 2 (28) 

<rr'> 

~ K} j COs(0i - 02 + 03 - 04 + 05 - 06)- 

o 

Here the second sum is over the hexagonal loops of the di- 
amond lattice (Fig. and the numbering inside the co- 
sine proceeds around the perimeter of the given hexagon. 
The diamond lattice is bipartite, so we define an orien- 
tation by declaring that links naturally point out of the 
"up-pointing" sites and into the "down-pointing" ones 
(corresponding to up- and down-pointing tetrahedra, re- 
spectively). We define an oriented electric field and vec- 
tor potential exactly as in Eqs. (|14I15|) . The Hamiltonian 



with the sense of the lattice curls determined as in the 
cubic case. Here (div6) r = 2n r , with n, = 0,±1. The 
Hamiltonian takes the form 

ftpd = ^^((curla)o- +vre° r ,) 2 - if ^cos(?r6 rr /). 
O* <rr') 

(33) 

One can derive an action in eigenstates of the dual vector 
potential as in the previous section, with the result: 

Z P d = E E ex P(-Spd), (34) 

and 

S pd = -Lhl(A) J2 (Ararr'+«rr-Qr'r) 
T,(rr') 

+ ^E(( curl «)o*+<-) 2 . (35) 

T,0* 

The constraint of magnetic charge quantization (div b) r £ 
2Z enters as before, giving rise to the temporal dual vec- 
tor potential fields in Z pc i- 
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III. EFFECTIVE THEORY 
A. Coulomb Phase Effective Action 

It is well-known that the low-energy description of the 
Coulomb phase of compact U(l) gauge theory is non- 
compact QED with no matter fields. In our case it will 
be convenient to work in dual variables to formulate the 
effective theory, which more naturally allows the inclu- 
sion of magnetic charge fluctuations. Depending on the 
purpose at hand, different formulations of the effective 
action will be useful. To fix notation, we show them all 
here. In each case the partition function is of the form: 

T r.fj, 

where a is a real- valued field. It will often be convenient 
to keep the full spatial lattice structure, in which case we 
can define a in terms of the microscopic variables by the 
temporal coarse-graining: 

Q v = W/ - — • ( 37 ) 

The brackets []/ denote an average over high-frequency 
modes in imaginary time, and we have subtracted the 
time-independent, non- fluctuating background a„i G 7rZ 
defined by: 

i(curla°) n . = -2e° r ,, (38) 

with aj? r = 0. If (as we always do in practice) we restrict 
attention to spatially periodic e® r , with zero average elec- 
tric flux in every direction, a can also be taken periodic. 
This subtraction simplifies the relation between the elec- 
tric field and dual vector potential, since 

ne rr , = ( curl(a - a°/2)) p . . (39) 

On a spacetime lattice, the action looks almost identi- 
cal to Scd in Eq. lf?7jl : 

S lat = Y E £(^W + ~ * 'r? (40) 
r (rr') 

T □* 

The only difference from the microscopic model is that 
the fields are now continuous. We will also have occasion 
to retain the spatial lattice structure but take the time- 
continuum limit, in which case we write 

S° tc = / drY^idr&n- + <5 rr - cv r ) 2 (41) 

(rr') 

+ ^ /drE(curl5) 2 . 



The parameters here are related to those in Eq. (|40ll by 
K, = (eg T ) -1 and U = g s /e, and <S rr now has units of in- 
verse time. This action corresponds to the effective dual 
Hamiltonian (in the sector with no magnetic charge): 

w° = ^E^ + |E( curi «) 2 ' ( 42 ) 

(rr') □» 

which can be obtained simply by expanding the cosine in 
Eq. I]22p. The lattice effective actions for the pyrochlore 
model look identical to Eqs. 1401 141(1 . where the only 
change necessary is the replacement of the sums over dual 
square plaquettes with sums over dual hexagons. 

Finally, in order to take the spatial continuum limit, we 
introduce a continuum 4- vector field T M (r, r), formally 
defined by the replacements 

& n -> Z(e; • T) (43) 
a rT — ► T T . 

Here are the vectors connecting the site r to its nearest 
neighbors, in either the cubic or diamond lattice, and I is 
a microscopic length on the order of the lattice spacing. 
Naively it seems we have thrown away too much informa- 
tion in the pyrochlore case, since there are four sites per 
unit cell, but only three spatial components of the con- 
tinuum vector field. However, there are also two gauge 
degrees of freedom per unit cell of the diamond lattice, 
corresponding to changing (diva) r on the two different 
sublattices. This leaves us with two transverse degrees of 
freedom, the same number as in the cubic model and the 
continuum theory. The full spacetime continuum theory 
is 

S °tc = 2^E/ dTd 3 r(d r ?i - d % T T f (44) 

i 

i<j J 

where IC C = KA and IA C —Ul (I is the microscopic length 
in Ea. 143(1 . The spatial/temporal components of T M have 
units of inverse length/time. The form is identical to the 
familiar (dual) Maxwell action for electromagnetism with 
photon velocity v p = \/U c K, c . This action can be used to 
obtain long-wavelength properties of the Coulomb phase 
of either microscopic model, but at the end of any calcu- 
lation the allowed spatial components of all vector fields 
are determined by the lattice structure. 

To calculate using any of these effective actions, it is 
convenient to implement a gauge- fixing procedure. The 
well-known manipulations of Faddeev and Popo^s^ tell 
us that we may add any function of the vector potential 
4-divergence to the action. The standard choice, in our 
nonstandard notation, is: 

Sfp = 2|^2 J drd 3 r(d T T T + v 2 p V ■ T) 2 . (45) 
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If we choose £ = 1, as usual the off-diagonal terms in 
the action are cancelled and we obtain the simple photon 
propagator: 



(T /1 (k,w n )T I/ (k', W ;)) = 



v 2 k 2 



(46) 

where the Euclidean metric is defined by g TT = ic, 
gri = 9ir = and gtj — Sij. While these expressions 
hold in the continuum, it is amusing to note that the 
obvious lattice regularization of Eq. I|45|l also cancels all 
off-diagonal terms on the cubic lattice. Working in the 
time-continuum action, we find the lattice propagator: 



(5 M (k,<x!„)a !/ (k',a;4)) 



-WiC ff (k) 



(47) 

where g(k) = 2^(1 — cosfci). On the diamond lattice, 
even more complicated Faddeev-Popov terms seem only 
to cancel some of the off-diagonal terms. 

Finally, it is important to note that this procedure is 
only legitimate for calculating expectation values of op- 
erators invariant under continuous dual gauge transfor- 
mations. Terms with only discrete dual gauge invariance 
(see Sec. IlllCl) must be handled more carefully; this issue 
arises in calculating the monopole propagator in Section 



B. RK Point Effective Action 

Given the microscopic gauge structure of our models, it 
is reasonable to conjecture that the low-energy effective 
degrees of freedom are simply noncompact U{\) gauge 
fields. The effective action should take the form of an 
expansion in the lowest-order terms involving the dual 
vector potential consistent with the symmetries. Gener- 
ically this would take the form of Eq. (|44ll : however, as 
our analysis of the microscopic models takes advantage 
of the special properties of the Rokhsar-Kivelson point to 
conclude that the "stiffness" or U c term for the electric 
field vanishes there, we must include a higher term. This 
leads us to propose that the RK point of both models is 
described by the effective action 

5 «* = 2)b ?/ drd3r( ^ T * - diT ^ )2 (48) 



^ J drd 3 r(Vx(VxT)) ; 



Here W c is the coefficient of the new term, which (when 
U c = 0) is marginal in the renormalization group sense 
and must be included. While the theory remains gauge 
invariant, the photon dispersion now vanishes quadrati- 
cally: ui ~ k 2 . The U c term is relevant, as can be seen 
by the usual power-counting procedure, so this action 
cannot describe a stable phase. A small positive U c will 
drive the system into the Coulomb phase, while a small 



negative U c will result in an electric field "crystal;" in 
the dimer language this is a state with long-range dimcr 
order and no flippable plaquettes. As we argue in Sec- 
tion precisely this picture describes the physics of the 
microscopic models near the RK point, with the relation 

U c OC (1 - V/Jring)- 

A very similar story is known to apply to the d = 2 
square lattice QDM 22 . In that case, electric-magnetic 
duality leads to a height model partition function, which 
describes the fluctuations of a discrete field living on the 
square lattice plaquettes in 2+1 dimensions. Just as 
above, softening the constraint of discreteness naturally 
leads to the Gaussian effective action: 



eight 



= / drd 2 r 



((d r /i) 2 + Ki(V/i) 2 + K 2 (V 2 /i) 2 ). (49) 



When ki = we have a description of the RK point, 
which has power-law dimer-dimer correlations. For small 
negative K\ the system goes into the staggered va- 
lence bond crystal state. These have natural analogs 
in the three-dimensional case. However, a small posi- 
tive K\ leads to a confining state with broken translation 



symmetr; 



.34.35 



via an instability of the Gaussian theory. 



This occurs because Eq. I|49[) with k\ > would describe 
the unstable Coulomb phase of pure 2+1-dimensional 
U(l) gauge theory. In three dimensions, however, the 
analogous phase is stable and should exist adjacent to 
the RK point. 

To calculate with Eq. I|48|l we again carry out a gauge- 
fixing procedure. We add the Faddeev-Popov term: 



Sfp-rk 



^3 J drd 3 r(a r T r +/C 2 V-T) 2 



(50) 



Although this does not cancel all off-diagonal terms, it 
leads to the relatively simple photon propagator: 



;T Al (k,cJ„)T i/ (k',^)) 

= (2ir) 4 6(Lo n +Lo' n )6(k + k')IC c M^(KL 



(51) 



where 
M TT 

Mij 
and 



M iT = 
1 

/(Wn,k) 
1 

/(w„,k) 



Klk 2 ) 



(52) 



[m 2 w 2 + kfk 2 + m 2 K, 2 c (k 2 - k 2 )} 
[hk,(k 2 - JC 2 c m 2 )] 



/K,k) = (uj 2 n +IC 2 c k 2 )(m 2 uj 2 n + k 4 ), 
with m = (K. C W C )- 1/2 . 



(53) 



Corrections 



While all gauge-invariant corrections to S® tc are irrel- 
evant in the renormalization group sense, in order to un- 
derstand the detailed realization of the U(l) spin liquid 
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(and the nearby RK point) in any microscopic model it 
will generally be important to consider corrections to the 
effective action. Furthermore, we need to specify rela- 
tions between the microscopic and effective degrees of 
freedom. Naively we could write ~ + oP^/2, but 
any corrections allowed by symmetry will generically be 
present. 

For ease of presentation, we focus on the cubic model; 
the lattice structure is unimportant for these results. The 
symmetries of the microscopic model are listed in detail 
in Appendix 1X1 It is important to note that, because the 
background a is not invariant under all lattice symme- 
tries, the dual vector potential transforms with additional 
shifts under these operations. Furthermore, the dual ef- 
fective degrees of freedom a and b transform exactly as 
their microscopic counterparts. 

Rather than attempt a painstaking enumeration of all 
allowed corrections to S° , we discuss some representative 
examples. The corrections naturally fall into two classes; 
those invariant under continuous dual gauge transforma- 
tions, and those not. In Hamiltonian language on the 
lattice, some typical terms in the first class are: 

H 1 = y ^(curlg) 2 + W^(curl6) 2 (54) 
□ □* 

<rr'> (rr') 

Terms involving discrete line integrals of both vector po- 
tentials around closed loops larger than single plaquettes 
are also allowed, as are terms containing the divergence 
of the fields. While these terms are irrelevant, they can 
have quantitative effects, presumably accessible in per- 
turbation theory. All such corrections are also irrelevant 
at the RK point, except the W term, which contributes 
to the W c term in Eq. (gHJ- 

More interesting are those terms lacking continuous 
dual gauge invariance. Here we only consider the single 
site terms; in the spacetime lattice action these take the 
form: 

oo 

S 2 — — 2_, E V2 i cos q (a° r ,) cos(2ga rr /) (55) 

T,(rr') q=l 
oo 

- ^2^2vl q cos{2qa rT ). 

T,r q— 1 

In the first term the a°-dependence is necessary to com- 
pensate the shifts in a under lattice symmetries. With 
appropriate choices of coefficients, when these corrections 
become large they have the effect of pinning the dual 
vector potential to take on discrete values. This allows 
us to interpolate explicitly between the effective theory 
Si at and the microscopic partiton function. Physically, 
the spatial V2 q part of S 2 is a magnetic charge hopping 
term. Its temporal companion is related to the discrete- 
ness of magnetic charge. Therefore, these corrections in- 
troduce magnetic charge fluctuations into the effective 



theory. As should be expected when magnetic monopoles 
are gapped, when these fluctuations are small there is no 
associated instability. Formally, this should be the case 
because the correlation functions of cos(2qa fl ) are local 
in space and time. For example, 

^ cos (2a rx (r)) cos (2a r , x (T r )) ^ cx 5„>5 TT >, (56) 

because continuous dual gauge invariance in S° makes 
each cosine into an independent random variable. 

Finally, we consider corrections to the scaling equali- 
ties between microscopic and effective degrees of freedom. 
Some representative contributions in the case of the dual 
vector potential are, in operator language: 

a„< ~ (<5v + oi„,/2) (57) 

oo 

+ E C2 i cos<1 ( a ° ' ) sm ( 2 9QV ) 

9=1 

OO 

+ E < m (^) 2n (curle) 2 D m+1 + --- 

n,m— 

This should really be interpreted as an instruction for 
what effective theory operators to use in the calculation 
of gauge invariant expectation values. These corrections 
are not important, and we will neglect them: the sec- 
ond term clearly only leads to local corrections, while 
the third term and others like it lead only to subdom- 
inant power laws, which will always be present on the 
lattice. The conclusion is the same for corrections to the 
magnetic field, but for completeness we show some of the 
representative terms: 

e ib -' ~ exp ( i b n , (58) 

OO 

+ * E < m (curlS) 2 n "(6 rr 2m+1 + ---) 

n,m— 

• j 1 + E cos 9 (a° r ' ) cos(2gq rr / ) + ■ ■ ■ j 

IV. PROPERTIES OF THE 17(1) SPIN LIQUID 

Using the effective theory developed above, we now 
treat some of the more striking properties of the U(l) 
spin liquid. Since the effective theories for both the 
generic Coulomb phase and the adjacent RK point are 
quadratic in the fields, all of the calculations can be done 
by simple Gaussian integrals. 

A. Excitations and Emergent Long-Range 
Interactions 

As already mentioned, the U(l) spin liquid supports a 
gapped S z = 1/2 spinon carrying electric gauge charge, 
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a gapped topological point defect that plays the role of 
a magnetic monopole, and a gapless, linearly dispersing 
"photon" with velocity v p . The spinon was discussed 
in Section Til CI and in the microscopic models takes the 
form of a tetrahedron (or octahedron) carrying nonzero 
S z . The gap to the spinons is very large, of order J z . 

The monopole is a classical configuration of the mag- 
netic field emanating from a point with nonzero magnetic 
charge at r . For a compact gauge field on, say, the cubic 
lattice, the classical configuration A of the vector poten- 
tial is given by minimizing 

E[A] = -^cos(7rB rr /) (59) 
("•') 

in the presence of a specified distribution of magnetic 
charge, which enters via the constraint (divZ?) r = 2<5 rro . 
In these expressions nB n > = (cm\A)u- Far from the cen- 
ter of the monopole the magnetic field will be small, the 
energy is approximately a sum over Bf r , , and the mini- 
mum is obtained by solving the (discrete) Laplace equa- 
tion of classical electrostatics. The compact nature of the 
theory only modifies the field near the core. 

In the effective theory, we can create a monopole with 
the operator 

m\ n = exp(iy^g rr /o: rr /). (60) 

(rr') 

The exact choice of B in this operator is somewhat arbi- 
trary; the requirement for monopole creation is that the 
surface integral of the magnetic field at large distances 
indicates the presence of one magnetic charge. This ar- 
bitrariness does not matter when calculating long-time 
monopole correlators, which will be dominated by the 
contribution from the lowest-energy monopole eigenstate 
(it is not hard to see that, in the effective theory, our 
operator has very good overlap with this state). As long 
as the magnetic field is chosen to spread the flux uni- 
formly from the monopole center, long-distance correla- 
tions should also be unaffected. Because magnetic charge 
is created at ro, this operator is not invariant under con- 
tinuous dual gauge transformations. The monopole gap 
is of order /C, and in the microscopic models is likely to 
be of order J r i ng as long as K, does not renormalize too 
much from its bare value. Therefore the monopoles have 
a much smaller gap than the spinons. 

Using the correspondence between microscopic and ef- 
fective variables, we can write down a creation operator 
that should have at least some overlap with the true mi- 
croscopic monopole eigenstate: 

m] o = exp(i ^2 A Tr >e Tr ,). (61) 

(rr') 

This is written in terms of the direct variables to connect 
with the spin language, where the monopoles are defect 
configurations in the XY component of the spin. Because 
of the half-integer electric field, the sign of rrv is changed 



by a 27r-shift of any one of the A rr ' ■ However, this overall 
sign does not fluctuate since the electric fields change 
only by integer steps. 

Both the monopoles and spinons feel an emergent 1/r 
interaction, even though the microscopic Hamiltonians 
contain only local operators. Consider the field due to 
a configuration of a few static electric gauge charges 
(the same discussion could be repeated for magnetic 
charge). This is given by solving Poisson's equation 
V-E = <Zi<5(r — r^) to obtain Coulomb's Law. By gauge 
invariance the longitudinal part of E does not fluctu- 
ate, so the field is simply given by its classical value plus 
transverse fluctuations. This zero point energy does not 
contribute to the energy difference between the charged 
state and the charge-free ground state, and one recovers 
the familiar 1/r Coulomb potential. In the RK point ef- 
fective theory, this interaction is absent for the spinons, 
but still present for monopoles, because there is no en- 
ergy cost for longitudinal electric fields. 

A word about the statistics of the charged excita- 
tions is in order. Certainly, both monopoles and elec- 
tric charges as discussed above are bosonic. However, a 
simple argument shows that a dyon, a bound state of elec- 
tric and magnetic charge, has Fermi statistics 3 ^. If such 
bound states exist and have lower energy than pure elec- 
tric charges, then one could say that there are fermionic 
spinons. While this is a nonuniversal question with no 
bearing on the low-energy physics, it can be relevant in 
understanding the possible transitions out of the U(l) 
spin liquid. The emergence of Fermi statistics in local 
bosonic models is also of great conceptual interest, al- 
though in the present case it is not clear what, if any, 
local microscopic interaction will bind electric and mag- 
netic charge. 

Finally we turn to the photon. In terms of the spins, 
it is a linearly dispersing S z — collective mode os- 
cillating between the Ising and XY parts of the spin 
vector. As for phonons, these gapless excitations make 
the contribution to the low-temperature specific heat 
C p hoton(T) cx (T/vp) 3 . Since in real materials it should 
be possible to quantitatively understand the phonon con- 
tribution to C(T) by measurement of elastic moduli 39 , 
this signature of the U(l) spin liquid should be easily 
accessible to experiments. Other potential, though likely 
more delicate, probes of the photon are low-temperature 
thermal conductivity and Raman scattering. Finally, the 
photon manifests itself in the power-law correlations dis- 
cussed in Section HV CI 



B. U(l) Topological Order 

In the theory of Z2-fractionalized phases, the notion 
of topological ground state degeneracy^ has been very 
useful both as a conceptual tool, and as a property that 
can be directly probed by experiments&i£*i£*i£. The de- 
generacy is associated with the topological sectors of a 
deconfining Zi gauge field, so it is natural to ask about 
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the generalization of these ideas to [/(l)-fractionalized 
states. 

We will work with the microscopic Hamiltonian of cu- 
bic frustrated gauge theory (Ea. 116(1 on a 3-torus, a cube 
with L lattice sites on a side and periodic boundary con- 
ditions: 



e [r+LeJj — e rj 

exp(ia[ r+Lei ] tj ) = exp(ia rj ). 



(62) 



We define an operator to measure the electric flux 
through a plane perpendicular to each of the three in- 
dependent cycles of the torus. For example, 



L-l 



(63) 



with similar definitions for These fluxes are con- 

stants of the motion, and are in fact conserved by any 
local, gauge invariant dynamics. Furthermore, because 
there is no electric charge in the ground state, by Gauss' 
Law we can make arbitrary incremental deformations of 
the specified plane without changing the flux. These 
fluxes define the electric topological sectors of the gauge 
theory^. 

Now imagine the model is in the Coulomb phase, where 
in the ground state $f = 0. We imagine threading 
one quantum $^ of electric flux through the system 
in, say, the ^-direction. Because the energy is propor- 
tional to J r E 2 , the most favorable situation is for the 
flux to spread itself uniformly through the system, so 
that E x — $|f /L 2 . This electric field is purely longi- 
tudinal and does not fluctuate, so there is a total en- 
ergy cost £ ~ U c (&q) 2 /L. Therefore, in the thermody- 
namic limit, all states with a finite number of electric 
flux quanta threading each direction of the 3-torus be- 
come degenerate with the ground state. Just as in the 
case of Z 2 topological order, these states are locally iden- 
tical to the ground state (local correlation functions will 
be unaffected), but globally distinct. 

What about magnetic topological sectors of the gauge 
field? Because [& rr /,$f] = 0, these can clearly be speci- 
fied simultaneously with the electric sectors. We define 
magnetic fluxes $f just as above: 



L-l 



(64) 



and similarly for &y z - This expression is a (lattice) sur- 
face integral over a plane bisecting the system, which has 
the topology of a 2-torus due to the periodic boundary 
conditions. Because this is a closed surface and b is the 
curl of a compact vector potential, the magnetic flux is 
quantized: $f = 2nf . 

However, there is an important difference from the case 
of electric flux: because magnetic charge fluctuates in 
the ground state, <&f is not a constant of the motion 




FIG. 6: Illustration of threading a magnetic flux quantum 
through one direction of a periodic system by creating a 
monopole-antimonopole pair (left), and moving it apart until 
the charges return to the same position and annihilate. One 
is left with no monpoles, and magnetic flux threading the 
system (right). 



in the microscopic theory. Within the Coulomb phase, 
though, the magnetic flux defined in the effective theory 
does commute with 7i°, the effective Hamiltonian writ- 
ten in dual variables. Magnetic charge fluctuations can 
mix states with different values of the flux by creating a 
monopole-antimonopole pair and separating the particles 
along one lattice direction, only to have them annihilate 
when they complete their periodic journey around the 
system (Fig. |SJ). Because the monopoles are gapped, the 
rate for this process is suppressed exponentially in the 
system size, and it does not occur for a large enough 
sample. As with the electric sectors, the magnetic flux 
sectors become degenerate ground states as L — > oo with 
energies scaling to zero as 1/L. In the microscopic the- 
ory, we expect that ($f ) will label the topological sectors 
within the Coulomb phase; in confining phases magnetic 
charge fluctuates wildly and $f is no longer a good quan- 
tum number, even approximately. 

To summarize, the U(l) spin liquid on a 3-torus has a 
topological degeneracy characterized by six integers, with 
energies that vanish as 



($f) 2 



£ 



U, 



+ K 



(4>o S ) 2 



(nf) 2 + «y + «) 



E\2 

y 



E\2 



(65) 



In the RK point effective action U c vanishes, so the elec- 
tric topological sectors have zero energy, with only cor- 
rections from irrelevant terms possible. 



C. Spin-Spin and Plaquette-Plaquette Correlators 

We first consider the two-point S z correlation function 
both in the Coulomb phase and at the RK point; the 
transverse part of the spin-spin correlator is not gauge- 
invariant and vanishes. The mappings of Section III CI 
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tell us this is given by the electric field correlator. In 
the cubic case we use the three-site unit cell containing 
{rir, ry, rz}. Because of the orientation convention for 
the electric field, we have — exp(iK • r)e r {, where 
Ko = (tt,tt,tt) and we take r = in the A sublattice. 
Therefore the spin correlation function will be shifted 
from the electric field one by Ko in the Brillouin zone. 
For the pyrochlore we use the unit cell of Section lll Al con- 
taining the four sites (R, i) , where R labels the centers of 
up-pointing tetrahedra (up-pointing diamond sites). In 
this case — eR^. 

We calculate the equal-time correlator C^(r — r') = 
(Ei (r)Ej (r')) in the gauge-fixed continuum theory for the 
Coulomb phase, where Ei — eijkdjTk- We illustrate the 
calculation with a single component: 



c£(R) 



K-c f d 3 \i fc 2 + ky ikR 
2v p J (2tt)3 k 6 



(66) 



= 7T~~ I 777~TT / d£l' sin 2 9' expfikR cos 7) , 
2v p J (27r) d J 

where we have imposed the hard momentum cutoff A, 7 
is the angle between k and R, 9 1 is the angle between 
k and the z-axis, and the O' integral is over the angular 
direction of k. Expanding the exponential in Legendre 
polynomials of cos 7, and using the addition theorem to 
rewrite these in terms of spherical harmonics, the integral 
can be done to find: 



Cf 2 (R) 



TT 2 V p R 4 



(2 cos 2 9 - 1). 



(67) 



Here 9 is the angle between R and the z-axis. We have 
dropped terms oscillating at the cutoff wavevector, which 
are unphysical artifacts of the hard cutoff. Alternatively, 
we can use a soft cutoff by inserting e~ fe / A in the in- 
tegrand, integrating over all fc-space, taking the limit 
R ^> A -1 and keeping only the dominant powers of R 
to recover the same result. In the general case we find: 



TT 2 V p R 6 



(2RiRj — R 2 Sij). 



(68) 



The striking angular dependence of this correlator is 
a manifestation of the inherent vectorial structure of the 
U(l) spin liquid, which comes in turn from the vector 
fields of U(l) gauge theory. For the pyrochlore model in 
the U(l) spin liquid phase we have: 

(SkiSk'j) ~ (eMej)iCi(R- R), (69) 

where are unit vectors connecting R with its nearest- 
neighbor diamond sites. In the cubic case: 

{Stflj) ~ exp (zK ■ (r - r'))Cf (r - r'). (70) 

In the RK point effective theory, a similar calculation 
shows that both the radial and angular dependence of the 
correlator changes. It takes the familiar "dipole" form: 



Cfj RK (R) - J (3RiRj - R 2 8 lj ), 



(71) 



falling off as 1/i? 3 . 

It is also interesting to calculate the correlations of the 
kinetic energy density. In our ring exchange models this 
naturally lives on the plaquettes and in the cubic case 
has the form 



1 



+ h.c.) = cos b„ 



(72) 



In the pyrochlore model Eq. (|72|l looks the same, but with 
ring exchange on the hexagonal plaquettes as in Eq. 
Long range order at k =/= in (e r i£ r >j) would indicate 
a plaquette density wave state with broken translation 
symmetry. To simplify notation we work out the cubic 
case, where 



(e ri e r >j) ~ (cos(6r»)cos(6r>j))o 

= exp(-(6 2 ) )cosh((6,A'j)o) 



(73) 



The prefactor involving (6 2 )o is a nonuniversal constant. 
At large separation, the second factor can be evaluated 
in the continuum theory, giving: 



M~Cl + (C«(r-r')) 2 + 



(74) 



We can immediately write down the magnetic field cor- 
relator by duality, which simply interchanges IA C and JC C 
within the Coulomb phase. Therefore 



cm) 



£ C^(R), 



(75) 



and the kinetic energy density exhibits power-law 1/i? 8 
correlations and a nontrivial angular dependence. 

The RK point theory does not have this self-duality 
property, so we need to evaluate another integral to find 
the magnetic field correlator. Using a hard cutoff, one 
finds only unphysical oscillatory terms. Using the e~ k ^ A 
soft cutoff, the result is Cff RK (R) oc 1/Ai? 6 . This 
cutoff-dependence indicates that the angular dependence 
is probably nonuniversal, since there should be other 
1/i? 6 contributions with nontrivial angular factors that 
have been left out of the continuum theory. 



D. Monopole Propagator 

In this section we calculate the monopole propagator 
and verify that it falls off exponentially in space and 
imaginary time, both in the Coulomb phase and at the 
RK point. This is one of the key properties of the U(l) 
spin liquid, since at the transition to any nearby confin- 
ing phase the monopolcs will proliferate and condense. 
It is most convenient focus on the equal-time monopole 
propagator (m R rriR'), working on the spacetime lattice. 
For R ^ R', this expectation value vanishes taken with 
respect to <S;° Qt because it creates magnetic charge at two 
points. To understand the true behavior we need to in- 
clude magnetic charge fluctuations and add the correc- 
tions Eq. I|55fl to the Gaussian action. 
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For simplicity, we consider only the restricted set of 
correction terms: 

—6S = V2 /J cos(2<S rT ) + v 2 2_, cos(af r /) cos(2a rr '). 

T,r r.(rr') 

(76) 

We want to calculate the propagator in perturbation the- 
ory in 5S: 



(4™R') 



(77) 



ex P(* B " ,5i "') eX P(- S lat - SS ) 



("') 

where B is the classical magnetic field due to a monopole- 
antimonopole pair at R and R', respectively. To obtain a 
nonvanishing contribution, we need to bring down correc- 
tion terms until the complex exponential has been mod- 
ified to create no magnetic charges. To lowest order we 
obtain: 

R' 

(m^TOR-) ~ 4 R ~\ exp(i^S rr -a rr ') JJ> cos(2a rr >)\ , 



<rr'> 



rr'=R 



(78) 

where the product of cosines is taken over the shortest 
path connecting R and R'; for simplicity we restrict our 
attention to geometries where this is unique, although 
this is inessential. We have also dropped the background 
field a , which gives only an overall sign. The gauge- 
invariant part of the expectation value comes from the 
term in the product of cosines that threads one magnetic 
flux quantum from R to R', giving a new magnetic field 
B' with zero divergence, and allowing us to express the 
result as the expectation value of a gauge-invariant op- 
erator: 



R-R'| 



(exp(^# r ,5 rr 0)o. (79) 

(rr') 



What the perturbation theory has done is exactly to con- 
nect the monopole and antimonopole with a Dirac string. 

The exponential decay is clear from the prefactor in 
Eq. I|79l) . To evaluate the corrections to this, we apply 
the Faddeev-Popov procedure to the functional integral 
for the gauge-invariant expectation value, and calculate 
using the photon propagator. We could not have done 
this at the outset since the original operator was not 
gauge- invariant. It is most convenient to integrate by 
parts inside the exponential to obtain a result in terms 
of the classical vector potential and the electric field: 

(exp(« Y - 4 rr' g rr')>0 = exp (-- ^ -KiA< 3 3 >0 



(rr') 



r,r' 



The necessary integral is quite difficult to evaluate, so we 
resort to power-counting to determine the largest pos- 
sible contribution. Consider a spatial separation of R 



between monopole and antimonopole. One contribution 
will come from the region near the pair, giving a factor 
of R e from the integration, R~ 2 from the two vector po- 
tential factors, and R~ A from the electric field correlator 
in the Coulomb phase. These factors multiply to give a 
constant, so we expect that the largest possible contribu- 
tion to the integral is logarithmic in R 7 which contributes 
only a power-law prefactor to the propagator. Other con- 
tributions involving regions far away from the pair make 
subdominant contributions. At the RK point there is an 
extra power of R from the electric field correlator, and the 
dominant possible contribution is linear in R and gives a 
correction to the correlation length. 

Finally it is clear from these considerations that the 
unequal-time propagator is not substantially different, 
and decays exponentially in space and time. 



V. EXACT GROUND STATE WAVEFUNCTION 

In this section we return to the microscopic spin mod- 
els, and use the exact ground state at the Rokhsar- 
Kivelson point to extract information about the physics 
nearby. We begin with a discussion of the structure of 
the wavefunction and some simple properties that can 
be seen analytically, then proceed to a discussion of the 
numerical evaluation of several quantities. In addition to 
various equal time correlation functions, we make use of 
a remarkable property of RK-type points discovered by 
Henleyai to approximately evaluate the imaginary-time 
monopole propagator. 

We will consider finite-size cubic and pyrochlore lat- 
tices in this section. To fix notation, for the cubic lattice 
L denotes the number of sites on a side of a periodic sys- 
tem. In the pyrochlore case we take periodic boundary 
conditons so that R and R + L&i are identified, and we 
can think of L as the number of tetrahedra along one 
periodic direction. 



A. General Properties 

The RK wavefunction is an equal-weight superposition 
(with positive coefficients) of all possible dimer cover- 
ings consistent with the local dimer constraint. In the 
spin models the electric field can only take on the values 
±1/2, and it is convenient to visualize dimer coverings as 
zero-divergence configurations of an electric field that can 
point only forward or backward on each link. Within the 
vast superposition there are components of every elec- 
tric topological sector. Since these are not mixed by 
the dynamics they should be thought of as degenerate 
but distinct ground states. The RK point exhibits the 
electric sectors of U(l) topological order, and as in the 
Gaussian RK action they are exactly degenerate even in 
a finite-size system. 

To understand what sector to focus on, it is profitable 
to consider first-order perturbation theory in SVNf away 
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from the RK point. The first-order shift in the ground 
state energy is 5E = 6V{Nf), where (Nf) is the average 
number of fiippable plaquettes in a given sector. For 
SV > 0, sectors with no fiippable plaquettes, and hence 
no dynamics, will have the lowest energy, and we reach 
the usual conclusion that a valence bond crystal obtains 
on this side of the RK point. For SV < 0, however, 
the sector with the greatest average flippability wins out. 
This is in some sense the most disordered sector, and 
should have <&f — 0; several numerical checks support 
this conclusion. We are interested in SV < 0, so we focus 
on the zero-flux sector. 

We also wish to consider sectors that violate the zero- 
divergence constraint and contain spinons. Within each 
such sector, the lowest energy state is again given by an 
equal-weight superposition of all electric field configura- 
tions consistent with a background charge density spec- 
ifying the spinon positions. These are exact eigenstates 
with energy J z /2 times the number of spinons, which in a 
finite-size system is always even. Consider in particular a 
sector with two spinons. There is clearly zero energy cost 
no matter how the spinons are moved around, so the RK 
point is deconfined and has no Coulomb potential be- 
tween static electric charges. Again, this is consistent 
with the effective action Eq. I|48|). 

Finally, we calculate the variational energy of the 
magnetic topological sectors. In general, we create a 
classical configuration of the vector potential by act- 
ing on the RK ground state with the operator O — 
exp(i J2(rr'} ^rr'e rr ')i j us ^ as f° r the special case of 
monopole creation in Eq. (|61() . Of course we expect only 
that the resulting state has some overlap with the desired 
eigenstate of the gauge theory. Denote the RK wavefunc- 
tion by 

where the sum is over all configurations of the electric 
field in the desired sector and Af is the number of states 
contributing to the sum. Defining Nf [{e}] to be the num- 
ber of fiippable plaquettes in a given configuration, and 
Pf[{e}, □] to be one if the specified plaquette is fiippable 
and zero otherwise, the variational energy of interest is 
given by: 

E var [A] = {iPrk\0 ] H rk O\<iP rk ) (82) 

-]L F /[{ e }' D ] ex P( 21 E Ar<e rr 

= E P/[{e},D][l-cos((curM) n )". 

{«,},□ 

Here we use notation appropriate to the cubic model; the 
generalization to the pyrochlore is obvious. 



Specifically, consider the case where O threads one 
quantum of magnetic flux through the x-direction in the 
cubic model. Then curl A — 2ir/L 2 for plaquettes in the 
x-direction and zero otherwise, and 

E var = n f L 3 J rmg (l - cos(27r/L 2 )) , (83) 

where nf is the average flippability per plaquette. For 
L — > oo this energy goes to zero as 1/L as expected, and 
further gives a rough value for the "magnetic stiffness" 
at the RK point: K, rs nfj r i ng . Using the numerical 
methods discussed below, for the cubic model in the zero- 
flux sector we find nf w 0.260. For the pyrochlore a 
similar calculation shows /C rs 2nfJ r i ng ; the factor of 
two arises because the flux passes through two kinds of 
hexagons. In this case we find numerically nf ~ 0.175. 

B. Monte Carlo Algorithms and Ergodicity 

Because the RK wavefunction has positive and equal 
weights, equal-time properties can be evaluated by infi- 
nite temperature Monte Carlo for the associated classical 
dimer model. The simplest possibile Monte Carlo step is 
to: 1) Randomly choose a plaquette. 2) If the plaque- 
tte is fiippable, flip it, otherwise do nothing. As desired, 
this algorithm preserves the electric flux $f. For the 
measurement of the quantities discussed below, the algo- 
rithm was run for as many as 10 11 Monte Carlo steps at 
a given system size. In some cases these long runs were 
necessary to achieve good accuracy because the desired 
quantity was very small. 

It is not clear whether the single ring move algorithm 
is ergodic within each electric field sector. For L = 2 cu- 
bic and pyrochlore systems we have performed an exact 
enumeration of all allowed configurations. In the cubic 
case the zero-flux sector contains 880 states; 864 of these 
have fiippable plaquettes (with n/ = 1/3 on average), 
and are connected under single ring moves. The other 16 
states have no fiippable plaquettes. For the pyrochlore 
model there are 384 zero-flux states connected under sin- 
gle ring moves (n/ = 1/4), and 12 zero-flux states with 
no fiippable hexagons. Evidently the single ring move 
algorithm is not ergodic within the zero-flux sector, but 
it does not matter if we fail to access states with no fiip- 
pable plaquettes. In our simulations we generate initial 
configurations for larger lattices by periodically repeat- 
ing a state from the fiippable part of the L — 2 zero-flux 
sector, so more insidious problems could occur if this sec- 
tor breaks up into multiple fiippable subsectors in larger 
lattices, each closed under single ring moves. 

For an analogous two-dimensional model, it is possi- 
ble to prove that single ring moves are ergodic in each 
electric flux sector (see Appendix [Bj . In the absence 
of a similar result in three dimensions, we have checked 
some of our results for the cubic model with an algo- 
rithm due to Barkema and Newman that we believe is 
probably ergodio^ 2 . This algorithm is illustrated in Fig- 
ure [7| in the present case the basic idea is to flip loops 
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FIG. 7: Illustration of the loop move algorithm of Barkema 
and Newman looking at one plane of the cubic model. One 
randomly chooses a starting point (A), and executes a ran- 
dom walk along the lattice bonds, moving only along the di- 
rection of the electric field vectors. Once the random walk 
intersects itself (B), the tail is removed and the electric field 
is reversed everywhere along the resulting loop, preserving 
the zero-divergence constraint. If only contractible loops are 
desired (as in our simulations within the zero-flux sector), 
one can simply throw away non-contractible loops and repeat 
these steps until a contractible loop results; this does not af- 
fect detailed balance. The single ring move only reverses the 
electric fields around elementary plaquettes like that at (C). 



of arbitrary length (including single rings), keeping only 
non-contractible loops that do not change <&f . All prop- 
erties measured using both algorithms gave the same re- 
sults; we present two examples in the following section. 
While we have not implemented an algorithm with "loop 
moves" for the pyrochlore model, due to the great simi- 
larity of all measured properties to their cubic analogs it 
seems unlikely that ergodicity is an issue. 

C. Equal Time Properties 

We measured the electric field-electric field correlator 
in both the cubic and pyrochlore models. To extract the 
asymptotic dependence most simply we focused on the 
correlators measured at a separation of half the system 
size. It should be noted that one cannot trivially extract 
the angular dependence discussed in Sect ion llV Cl bv this 
method, since the two electric field vectors in the correla- 
tor will be connected by many paths with length scaling 
as L but with different angles. For the cubic model we 
show our results for (e x (L~x./2)e x (Q)) in Figure |S| Data 
from both algorithms are shown and the agreement with 
the 1/R 3 decay predicted by the effective theory is very 
good. The same conclusion obtains for the pyrochlore 
model (but with only the single ring move algorithm) . In 
that case the correlator measured was (eo(£ai/2)eo(0)); 
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FIG. 8: Log-log plot of the equal time electric field correlator 
at the RK point of the cubic model, in the orientation dis- 
cussed in the text. The circles denote data from the single 
ring move algorithm, while the diamonds were obtained by 
the loop move algorithm. Error bars are on the order of the 
symbol size. The line is a guide to the eye with slope —3 to 
show the very good agreement with the 1/R 3 decay expected 
from the effective action. 



the data are shown in Figure Measurements of other 
orientations of the electric field correlator in both models 
all showed the same 1/R 3 decay. 

To measure the equal time monopole propagator it 
is necessary to address some subtleties that arise when 
putting objects with gauge charge in a finite-size sys- 
tem. First, it is clear that because the system is a com- 
pact manifold with no boundary, it must have zero total 
magnetic charge. Since the operator mJmR/ creates a 
monopole-antimonopole pair, this is not a problem for 
the propagator. However, consider the cubic lattice and 
let R — R' = n x x ^ 0. For any value of n x there will 
be planes where & x ^ ^irn^ , which is not allowed for a 
compact vector potential as discussed in Section IIV Bl 
We can ameliorate this problem by creating a double 
strength monopole-antimonopole pair, but even in this 
case the only allowed separation is n x — L/2. Because of 
these complications, in the cubic case we only measure 
the double-strength propagator at a separation of half 
the system size: 

Gf(L/2)^((m; R+ix/2) ) 2 (m R ) 2 ). (84) 

By time-reversal or Ising symmetry, G^f is real. We ob- 
tained the vector potential appropriate for Eq. (|84|l by 
first solving a discrete Poisson's equation by numerical 
matrix inversion, then feeding the solution into Eq. I|59|) . 
This was minimized by a combination of simulated an- 
nealing and direct minimization (i.e. zero-temperature 
simulated annealing). Because of the very rapid decay, it 
was only practical to measure G^f for L < 8; results are 
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FIG. 9: Log-log plot of the pyrochlore RK point electric field 
correlator in the orientation discussed in the text. As in Fig. 
|5| the line is a guide to the eye with slope —3. 



shown in Figure I1UI and are consistent with exponential 
decay with power-law corrections. 

On the pyrochlore lattice it is not possible to have 
double-strength magnetic charges in the microscopic 
model (see Sec. Ill E|) . so for consistency with flux quan- 
tization we must consider a more complicated geometry 
We look at the following propagator, which creates two 
monopoles on adjacent up- and down-pointing (dual) di- 
amond sites, and similarly two antimonopoles separated 
by a distance L/2: 



G™(L/2) 



(85) 



, '(R+La /2) m (R'+La /2) ? 



where R is an up-pointing diamond site and R' is the 
down-pointing site directly above. The results are shown 
in Figure ITT1 and are consistent with exponential decay; 
we do not believe the apparent lack of substantial power- 
law corrections is significant, since the coefficient of these 
corrections is presumably nonuniversal. 

In the cubic model, we also measured the potential 
between a pair of static spinons in first-order perturba- 
tion theory away from the RK point. In the notation of 
Sec. IV Al we consider SV < 0. The coefficient of \8V\ in 
the spinon-spinon potential is given by: 

V S pinon(R-) = (Nf (no spinons) - N f (spinons)). (86) 

To measure this quantity, we generated sectors with 
one S* = 1/2 spinon at r = 0, and another at r = 
(L/2)x+ (L/2 — l)y for L > 4. The results are shown in 
Figure lT2l and are consistent with al/r Coulomb poten- 
tial. While the perturbation theory does not necessarily 
"know" whether a deconfining phase with a Coulomb po- 
tential can be stable, it does indicate U c > for small 
SV < 0. Since the Coulomb phase is stable, we conclude 
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FIG. 10: Cubic model equal-time monopole propagator. The 
vertical axis is the logarithm of the absolute value of the 
double-strength propagator discussed in the text, and the hor- 
izontal axis is system size. Data for both the single ring move 
(circles) and loop move (diamonds) algorithms are shown. For 
the loop move data the error bars are smaller than the symbol 
size, and for the first two single ring move data points they 
are obscured by the other symbols. The results are consistent 
with exponential decay with power-law corrections. 
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FIG. 11: Pyrochlore equal-time monopole propagator. The 
vertical axis is the logarithm of the absolute value of the prop- 
agator Eq. I|85|l . The results are consistent with exponential 
decay, with no indication of power-law corrections. 



it indeed exists adjacent to the RK point over a finite 
range > SV > SV min . 



D. Monopole Gap 

In our models it is possible to approximately measure 
quantum imaginary time correlation functions using only 
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FIG. 12: Plot of the spinon-spinon potential defined in 
Eq. 1861 for several system sizes of the cubic model. The 
horizontal axis is the distance between the spinons, R = 
yi 2 /2-L + l. The curve is a fit of the data to the func- 
tional form Vi + V 2 /R. 



the classical Monte Carlo dynamics. The key observa- 
tion, due to Henley^i, is that at the RK point the master 
equation describing the Monte Carlo dynamics of the sin- 
gle ring move algorithm is identical to the imaginary-time 
Schrodinger equation. We will be interested in monopole 
correlation functions built from the equal-time correla- 
tors constructed above. In the cubic case, for example, 
we measure: 



C c M (r) = (0 M exp(- 



-tHrk)O m 



where 0\ { — 



= 



jt 

l (R+Lx/2) 



2/ N 2 



(87) 

The analogous defi- 



nition of C^ 1 (t) for the pyrochlore model is constructed 
from the equal-time propagator Eq. 

To extract information about the spectrum from the 
simulation, it is useful to relate the units of classical 
Monte Carlo time to quantum imaginary time. Let N p 
be the total number of plaquettes in the system. If we 
begin with some electric field configuration {e} and exe- 
cute one Monte Carlo step, the probability of remaining 
in the same state is (1 — Nf[{e}]/N p ), while that of mak- 
ing a transition to each of the states accessible by flip- 
ping one plaquette is 1/N P . We now consider the time- 
discretized imaginary-time Schrodinger equation and fix 
the timestep At to recover the same values, which now 
enter as probabillity amplitudes. We consider: 

exp(-ArH RK )\{e}) (88) 

= (l- ATH R K + 0{(ATH RK ) 2 ))\{e}) 

= (l - J ri n fl AriV / )|{e}) + J rmg ATj2 \W}), 

where the final sum is over those electric field configura- 
tions connected to {e} by a single ring move. It is only 
valid to neglect the higher-order terms in AtTLrk when 
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FIG. 13: Logarithmic plot of the approximate monopole- 
antimonopole imaginary time correlator C c (r) for the cubic 
model. The horizontal axis is imaginary time in units of J~ ng ■ 
Data is shown for L — 4,6, 8; as discussed in the text, as L 
increases the exponential decay in imaginary time becomes 
cleaner, suggesting the monopole is indeed gapped. 



Nf/Np <C 1; in fact, the correspondence between the 
discrete-time quantum dynamics and the Monte Carlo 
dynamics (which are necessarily discrete in our simula- 
tion), is only strictly valid in this limit. In our models 
Nf/Np is of order 1/5, so we expect to at best make 
quantitative errors of about 10%, and at worst to get the 
wrong answer. This problem cannot be alleviated by the 
means at hand because the discrete classical algorithm 
fixes the time-step for the quantum dynamics. A more 
careful treatment would require a direct simulation of 
the master equation with control over At, or quantum 
Monte Carlo. 

If we nevertheless expand the exponential, setting 
At = (JringNp) -1 correctly matches the classical prob- 
abilities and quantum amplitudes. Using this relation, 
we write down the appropriate quantity to simulate to 
measure the imaginary time correlator: 



C (t 1 1 JringNp) 



{e} 



1 



N 



T(t,{e}) 



(89) 



£ (G M ({e}))*0 M (T(t,{e})) 

T(t,{e}) 

Here t is an integer number of steps in Monte Carlo time, 
and T(t, {e}) labels all possible Monte Carlo time evolu- 
tions of time t starting from the configuration {e}. The 



number of such evolutions is N 



T(t,{e}), 



and Af is the total 



number of electric field configurations in the zero-flux sec- 
tor. Finally, we have defined 6 M \{e}) = C M ({e})|{e}). 

We have measured the monopole imaginary time cor- 
relators for both the cubic and pyrochlore models; re- 
sults are shown in Figures ED an d \U\ Note that the 
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FIG. 14: Logarithmic plot of Cp'(r), the approximate 
monopole-antimonopole imaginary time correlator in the py- 
rochlore model. The horizontal axis is imaginary time in units 

°f j 'ring • 



extreme long time behavior is not related to the gap; in- 
stead, O m \iPrk) will have some overlap with the ground 
state, causing C A/ (r) to approach a constant as r — > oo. 
As we go to larger system sizes and the spatial separa- 
tion of the monopole-antimonopole pair in the correlator 
increases, there is a clean exponential decay persisting 
for longer times. This is illustrated for several system 
sizes of the cubic model in Fig. 1131 From the numerical 
data one can extract a rough value for the energy of the 
monopole-antimonopole configuration put in at t = 0; 
in both cases this is about A.hJ r i ng . While these results 
are approximate, it is unlikely that the monopoles are 
somehow gapless given the circumstantial evidence for a 
gap, particularly the very rapid exponential decay of the 
equal-time propagator. 



VI. DISCUSSION 

In this paper we have argued for the existence of a 
novel fractionalized quantum disordered state, the £7(1) 
spin liquid, in spin models with a global £7(1) symmetry. 
While we found it necessary to add an extra interaction to 
the easy-axis limit of the pyrochlore Heisenberg model to 
proceed analytically, we speculate that this may be only 
a crutch and that the pure easy-axis model is in the U(l) 
spin liquid phase. Both the pyrochlore and cubic models 
have no sign problem, so this issue could be resolved by 
quantum Monte Carlo simulations 43 . While it would be a 
remarkable result if the easy-axis pyrochlore Heisenberg 
antifcrromagnet were a £7(1) spin liquid, we believe the 
results in this paper are enough of an "existence proof" 
that the phase is likely to be present in less tractable but 
more realistic microscopic models. 

The converse problem, of determining the general re- 
quirements for a spin model to exhibit an emergent £7(1) 



gauge structure, is of course more challenging. It would 
appear that, since monopole excitations are required to 
enable a transition to a neighboring confined phase, any 
such gauge theory must be compact, and hence appear 
on the lattice. Thus, unlike emergent global symmetries, 
which can appear at low energies and long wavelengths, 
an emergent £7(1) gauge structure would appear to admit 
only a limited degree of spatial coarse-graining. Given 
the additional need for spin-carrying degrees of freedom, 
the unit cell of an emergent gauge theory in a pure spin 
model would appear to necessarily contain more than 
a single S — 1/2 spin. Based on the present examples, 
which take advantage of the bipartite structure of the cu- 
bic and diamond lattices, it is tempting to speculate that 
the models should contain a natural bipartite sublattice. 
It is easy to see, however, that the emergent gauge struc- 
ture in these models is stable to adding arbitrary weak 
additional multi-spin interactions, including those that 
break any bipartite symmetries. Further understanding, 
including the important issue of whether the £7(1) spin 
liquid is possible in a magnet with global SU(2) sym- 
metry, and if yes, of what its properties might be, must 
await further study. 

Certainly, the most interesting issue is whether the 
£7(1) spin liquid can be found in a real material. The 
simplest signature is the additive T 3 contribution to the 
specific heat from the photon. Because it is likely to be 
possible to quantitatively understand the phonon specific 
heat via independent measurements^, this could provide 
a relatively clear and simple test for emergent photons. 
Further theoretical work may be necessary to understand 
more delicate probes of the photon, like heat current and 
(possibly) Raman scattering; phonon- "photon" interac- 
tions and disorder may play a nontrivial role. While dis- 
order may be important for some experimental proper- 
ties, we remark that simple arguments demonstrate that 
disorder does not destabilize the £7(1) gauge structure or 
the gapless photon. 

Since topological order can coexist with conventional 
broken symmetry, spin liquids are not the only good 
places to look for emergent photons. One can imag- 
ine, for example, condensing spin-carrying but gauge- 
neutral excitations; in our models this would lead to an 
ordered magnetic state with gapless spin waves and a 
£7(1) gauge structure. Understanding the possibilities 
for and properties of phases near the £7(1) spin liquid is 
another problem worthy of consideration if contact is to 
be made with experiment. Very recently, Senthil, Vo- 
jta and Sachdev have made the interesting suggestion 
that £7(l)-fractionalization may provide the answer to 
some of the puzzles of heavy fermion materials^. We 
leave these intriguing issues aside to remark that £7(1)- 
fractionalized states are a remarkable and novel possibil- 
ity for the physics of strongly correlated electrons in three 
dimensions, heretofore a relatively unexplored area. 

Finally, we note that Huse, Krauth, Moessner and 
Sondhi are considering a model closely related to ours 45 . 
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APPENDIX A: ENUMERATION OF 
SYMMETRIES 

We here enumerate the action of the various symme- 
tries of the cubic model on all the microscopic operators, 
both in dual and original variables. 

• Discrete translations: r — > r + R. The electric field 
transforms as e r i — > er+R^, with similar expres- 
sions for b and a. Because of the non-translation- 
invariant background a , the dual vector potential 
obeys the more complicated transformation law: 

a xi — » a R+r j + 5a R+r ^, (Al) 

where the shift 5a is defined to satisfy 

±-(curl8a) r+Rii = e r+Rti -e° ri . (A2) 

• Lattice rotations: Let TZ rotate the lattice into it- 
self about some fixed origin. For example, take the 
origin to be r = and make a 7r/2-rotation about 
the z-axis. Then lZ{n x x + n y y + n z z,x] = [— %x + 
n x y + n z z, y\, and so on. We have e r i — > e-ji[ r i], an d 
similarly for b and a. Again a X i — ► a-ji[ r i\ + 5a-R[ r i] > 
with (cml5a) n[ri] = n(e^ [ri] - e°J. 

• Reflections: Let T denote reflection about a plane 
with normal vector x, y or z. For definiteness, 
choose the plane with normal x at x — 0, then: 

e rx — » er(r),-x = -e^( r )-x,x (A3) 

for i = y,z. The vector potential obeys a similar 
transformation law. The magnetic field and dual 
vector potential arc pseudovectors and transform 
under reflections with an additional minus sign. 
Again, a transforms with an appropriate shift 5a 
to compensate for changes in the background under 
reflections. 

• Gauge and dual gauge invariance: Under gauge 
transformations, a rr > — > a rr > + Xr' — Xr, where x is a 



phase defined on the cubic sites. Similarly we have 
the dual gauge transformations a n > — > a n >+\ r r — A r , 
with A r € 7rZ. In the action the dual gauge trans- 
formations can be spacetime-dependent, with 

a rl {r) -> a n (r) + A r+ei (r) - A r (r) (A4) 
a XT -> a rr + A r (r + e) - A r (r). 

• Ising (or particle-hole) symmetry: e — > — e, a — ► —a 
and 6 — > —6. Also a — > a — a. 

• Time reversal: Since spin operators obey S — > — S 
under time-reversal, we have: e — > — e, a — > a + 7r 
and 6 — > 6. For consistency with the electric field 
transformation law, a — > a — a. In the action the 
latter relation continues to hold, as long as we also 
send r — > — r. Furthermore, the temporal com- 
ponent of the dual vector potential transforms as 
a rr (r) — > — a rr (— t — e). Note that the situation 
here is the reverse of that in real electromagnetism, 
where B changes sign under time reversal and E is 
invariant. 

These symmetries are all present in analogous forms in 
the diamond lattice gauge theory. There is also an addi- 
tional global symmetry, which exchanges the two sublat- 
tices of up- and down-pointing sites. 

APPENDIX B: ERGODICITY OF SINGLE RING 
MOVES IN A SQUARE LATTICE MODEL 

We consider the classical dimcr model on the square 
lattice at infinite temperature, with two dimers touching 
every site. This is the d = 2 analog of the RK point 
of our cubic model. We work in the electric field lan- 
guage, where e r i = ±1/2 with i = x,y, and (dive) r = 0. 
Consider an L x L system (L even) with periodic bound- 
ary conditions: e\ T +Lei],j = e rjj where e, = x, y. The 
zero-flux sector is specified by the condition: 

L-l L-l 

y ' e n x x,j/ = y ' e n y y,x = 0. (Bl) 
n x — n y —0 

We will show that any two states in this sector are con- 
nected by a sequence of ring exchange moves on single 
square plaquettes. 

The key step is to go to a height representation on 
the dual lattice with sites at the plaquette centers r = 
r + (x + y)/2. We define: 

ft(r + x)-ft(r) = -2e [r+(x _ y)/2 ] 1 j / (B2) 
h(r + y)-h(r) = 2e [r+(y _ x)/2];a: 

The content of this definition is that the height in- 
creases/decreases by one if we cross a link of the direct 
lattice with electric field pointing to the right/left. If wc 
fix the value of the height on one site of the dual lattice 
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these definitions determine it uniquely everywhere. The 
height is well-defined because the electric field has zero 
divergence, and can consistently be taken to have peri- 
odic boundary conditions because we are in the zero-flux 
sector. 

A flippable plaquette has a height either above or be- 
low all of its four neighbors, depending on its orientation. 
Therefore, in a given configuration, the plaquettes with 
minimum and maximum height will always be flippable, 
and every configuration has at least two flippable pla- 
quettes. There are also two states with every plaquette 
flippable; up to overall shifts of the height these have 

/i(n a x + % y + (x + y)/2) = ±|(l + (-!)"-+"»). (B3) 



These two configurations are connected by single ring 
moves, since we can flip all the maximum height plaque- 
ttes in one to go to the other. 

To complete the proof, we label configurations by 
Ah = h max — h rn i n . This clearly takes on the minimum 
possible value of unity in only the two maximally flip- 
pable states. Suppose we are in some other state with 
Ah > 1. Then we can flip plaquettes of maximum height 
until Ah is reduced by 1. This procedure can be repeated 
until Ah = 1 and we reach one of the maximally flippable 
states. We have thus shown that any two states are con- 
nected by a sequence of single ring moves, because the 
maximally flippable states are connected to all states and 
to each other. 
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